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We consider excitations with the quantum numbers of a hole in the one dimensional Hubbard 
model below half-filling. We calculate the finite-size corrections to the energy. The results are then 
used to determine threshold singularities in the single-particle Green's function for commensurate 
fillings. We present the analogous results for the Yang-Gaudin model (electron gas with ^-function 
interactions). 
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I. INTRODUCTION 

The Hubbard model constitutes a key paradigm for strong correlation effects in one dimensional (ID) electron 
systems 1 . Its Hamiltonian is 
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where n,j t<T — ct ^Cj-.o- and rij = n^-j- + In the following discussion the magnetic field B will be set to zero, but we 
will reinstate it in the calculations in sections III A-V. The Hubbard model is solvable by Bethe Ansatz 2 and many 
exact results are available in the literature 1 . Of particular interest in view of experimental applications are dynamical 
response functions such as the single-particle spectral function 
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The spectral function is measured in angle-resolved photoemission experiments. Such measurements on the quasi-lD 
organic conductor TTF-TCNQ have been interprested in terms of A(u,q) of the ID Hubbard model 3 ' 4 . While high 
quality numerical results are available from dynamical density matrix renormalization group computations 4,5 , it is 
so far not possible to calculate (2) analytically from the exact solution. However, using a field theory approach it 
is possible to determine low-energy properties exactly. In particular, the singularity as a function of u> at the Fermi 
wave number can be obtained using Luttinger liquid theory 6 . The low-energy physics of the Hubbard model in zero 
magnetic field is described by a spinful Luttinger liquid with Hamiltonian H = H c + H s , where 1,7 
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irrelevant operators. 
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Here K s — 1 (we are concerned with the B = case for the time being) and the spin and charge velocities v c s as well 
as the Luttinger parameter K c are known functions of the density and interaction strength 1 . The Bose fields $ a and 
the dual fields Q a fulfil the commutation relations 
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The spectrum of low-lying excitations (relative to the ground state) in a large but finite system of size L is given 
by 1,23 
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where AiV Q , D a and N± are integer or half-odd integer "quantum numbers" subject to the selection rules 
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We note that the form of the finite-size corrections implies that in the finite volume the spin and charge sectors are not 
independent but are in fact coupled through the boundary conditions of the fields <£> Q , 9 Q . Neglecting the effects of 
the irrelevant operators in (3) makes it possible to calculate A(u>, q) at low energies 1,6,7 . The spectral function is found 
to exhibit singularities following the dispersions of the collective spin ("spinon") and charge ("holon") excitations. 
The exponents characterizing these singularities are given in terms of the quantum numbers A7V Q , D a and 7V^ . 




FIG. 1: (a) Support of the single-particle spectral function. A(u,P) is non-vanishing in the shaded areas. is the 

dispersion relation of the collective spin excitations ("spinons"). The blue lines depict the linearized dispersions ±v s (P — fcj?) 
and ±v s (P + /of) that underlie the Luttinger liquid approximation, (b) Structure of the single-particle spectral function in 
the Luttinger liquid approximation for momenta close to kp and P < kp. There is a singularity at negative frequencies 
uj = —Vs(kF — P) and a power-law "shoulder" at positive frequencies lj = v c (kF — P). 



In a series of recent works 8-19 it was demonstrated for the case of a spinless fermions, that neglecting the irrelevant 
operators perturbing the Luttinger liquid Hamiltonian leads in general to incorrect results for singularities in response 
functions. Using a mapping to a Luttinger liquid coupled to a mobile impurity and taking the leading irrelevant 
operators into account non-perturbativcly it is possible to determine the exact singularities in response functions 8-20 . 
Crucially, these singularities are generally momentum dependent. Two recent preprints have addressed the general- 
ization to spinful fermions 21,22 . In particular, Ref. 22 derives expressions for the exponents /i n ,± characterizing the 
the singularities of the single-particle spectral function (2). The resulting spectral function is depicted in Fig. 2. 




FIG. 2: Power-law singularity in the single-hole spectral function for \P\ < Hf- There is a singularity at negative frequencies 
uj — e s (P) with a momentum-dependent exponent /io,- and a power-law "shoulder" at positive freuencies to = — e s (P). 

Our goal is to provide an exact expression for the threshold exponents /io,- in terms of the microscopic parameters 
entering the Hamiltonian, i.e. U, [i and B in the case of the Hubbard model. This is achieved by determining 
the finite-size corrections to low-lying energy levels in presence of a high-energy excitation. Comparing the results 
obtained from the Bethe Ansatz solution of the Hubbard model as well as the closely related Yang-Gaudin model to 
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field theory predictions, we are able to derive explicit expressions for the threshold exponents. In the Yang-Gaudin 
case our explicit results agree with the relations of spectrum and exponents proposed in Ref. 22 for Galilei- invariant 
models 31 . 

There is one caveat for the case of the Hubbard model. As for any lattice model it is possible to generate low-lying 
excitations for any momentum by combining multiple Umklapp processes if the band filling is incommensurate. A 
completely analogous situation is encountered for spinless fcrmions 18 . Hence, for incommensurate fillings no thresholds 
exist. While we would still expect the spectral function to feature peaks asociated with the thresholds of particular 
excitations, in particular those involving small numbers of holons and spinons, these peaks will no longer correspond 
to singularities. In order to circumvent this problem we will consider only the case of commensurate fillings in the 
Hubbard model. 

The outline of this paper is as follows. In section II we briefly review the Bethe Ansatz description of the ground 
state of the Hubbard model. In section III we present the excitations that give rise to low-energy thresholds around 
the Fermi momentum in the Hubbard model. In sections IV and V we determine the finite-size spectrum of excited 
states describing these thresholds. In section VI we relate these results to the field theory treatment of the threshold 
problem and extract threshold exponents. In section VII we summarize the analogous results for the Yang-Gaudin 
model. 

II. BETHE ANSATZ EQUATIONS FOR THE HUBBARD MODEL 

The logarithmic form of the Bethe ansatz equations of the Hubbard model for TV electrons out of which M have 
spin down (for real solutions only) is 

k j L = 2.i j -^e( s ^zJ±\ j = 1 ,..., N , 

h{ bs =^)-^t»{^ £ ) .°=V..,M. (7) 
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the length of the lattice L is taken to be even, 9(x) = 2arctan(a;) and Ij, J a are integer or half-odd integer numbers 
that arise due to the multivaluedness of the logarithm. They are subject to the "selection rules" 

j . f integer if M is even 

J | half — odd integer if M is odd, ^ ' 
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if N - M is odd 
odd integer if N ~ M is even, 
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The energy (in units of t) and momentum of such Bcthc ansatz states are 
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A. Ground State below half-filling 

Following section 7.7 of [1] we now take consider the ground state below half filling. We have 

N = N GS , M = M GS , (13) 
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where we take Nqs = 2x odd and Mqs = odd. We note that in zero magnetic field we have Mqs = ^§ s -- Our choice 
for N and M implies that Ij are half-odd integers and J a are integers. In the ground state all vacancies are filled 
symmetrically around zero 

Ij =3 2 ' 3 = ' ' ' ' ' GS ' 

M GS + 1 

J a = a ^| , a = 1, . . . , Mqs ■ (14) 

It follows from (12) that the ground state momentum is zero. The bulk ground state energy can be expressed in terms 
of the solution of the following set of coupled integral equations 

1 f A 

p c (k) = h cosfc / dA ai(sin/s - A) p s (A) , (15) 

2tt J-a 
rQ rA 
p s (A) = / dk ai(A-sinfc) p c {k) - / dA' a 2 (A - A')p s {A') , (16) 
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The integrated densities yield the total number of electrons per site and the number of electrons with spin down per 
site respectively 

dk Pc (k) = ^=n c , [ A dAp s (A) = M?* = ^± = n s . (18) 

i L J-A ^ 

The integration boundaries Q and .4 can be fixed in terms of A^gs and Mqs hy these equations. Alternatively one 
can define dressed energies by 



e c (k) = -2 cosfc - p - 2u - B + / dA ai(sinfc — A) e s (A) , 

J-A 

rQ fA 

e s (A) = 2B+ dk cos(fc) ai(sin/c- A)e c (fc) - / dA' a 2 (A - A') e s (A') . (19) 

J-Q J-A 

Here the integration boundaries ±Q and ±A are by definition the points at which the dressed energies switch sign, 
so that they are determined as functions of the chemical potential and the magnetic fields via the conditions 

e c (±Q) = 0, e s (±A)=0. (20) 

The bulk ground state energy per site is 

/Q rQ (ifo 

dk (-2 cosfc - p - B - 2u) p c {k) + 2Bn s + u = / — e c (fc) +u. (21) 
-O J-Q 27r 



III. EXCITATIONS WITH CHARGE -e AND SPIN \ 

We now consider an excitation over the ground state with the quantum numbers of a hole with spin down. As the 
total charge must be one less than in the ground state we must have 

N = Nqs ~ I- (22) 

Recalling that the z-component of total spin quantum number is S z = N ~ 2M we see that relative to the ground state 
we must have 

2/ 4- 1 

5S Z = — ^— , I G No- (23) 

Equation (23) requires an explanation. As shown in 24 the Bethe ansatz states only provide highest-weight states of 
the SO (4) symmetry 25 of the Hubbard model. The corresponding SO (4) multiplet is then obtained by acting with 
lowering operators. As a result, for any excitation given by the Bethe ansatz with 8S Z = =^ for non-negative integers 
I we can construct an excited state with 5S Z = \ by acting with the spin lowering operator. In presence of a magnetic 
field this shifts the energy by —2BI. 
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A. Holon-Spinon Excitation 

The simplest excitation with the quantum numbers (23), (22) is two-parametric 1 and obtained by setting 

iV = JV GS -l, M = M GS -1, (24) 

in the Bethe ansatz equations (7). It then follows from (9), (10) that both Ij and J a are integers. In order to see 
that we are dealing with a two-parametric excitation we consider the number of vacancies for the integers Ij and J a 
(11). As 



L L 
2 J ~ 2 



(25) 



there is the same number of vacancies for the Ij's as in the ground state, but we have one fewer integer, i.e. one 
additional "hole" . Similarly we have 



\Ja\ < ^(N GS -M GS -1), 



(26) 



which tells us that we have the same number of vacancies as in the ground state but one root less, which leaves one 
"hole" . It is useful to plot the distribution of integers. This is done in Fig. 3. The energy and momentum of the 
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FIG. 3: Distribution of integers for the holon-spinon excitation for Ags = 14, Mgs = 7. The positions of the two holes are 
I h and A h respectively. We have shown the configuration where the Ij' range from -^rp — 1 to -^rp, but we equally well can 
choose them -^rp < Ij < — £p — 1. The two possibilities give rise to the two signs in (27). 



"holon-spinon" excitation described above are determined in Chapter 7 of Ref.l 

E hs = -e c (k h ) - e s (A h ) , 

Phs = -p c (k h )-p s (A h )±Tm c , 

where the dressed energies are defined in (19) and the dressed momenta are 



(27) 



p c (A) - 2tt / dk' p c {k') , 



p s (A) = n(n c - n s ) - 2ir / dA p s (A). 

Ja 



(28) 



The extra contribution ±irn c in (27) arises as the distribution of i^'s is asymmetric with respect to the origin. In zero 
magnetic field the equation for p s (A) can be simplified 



Ps(A) 



B=0 



dk arctan 



exp 



(-|-(A-sin(£0)]p c (fc) 



(29) 



In Fig. 4 we plot the upper and lower boundaries of the holon-spinon continuum for two different band fillings in zero 
magnetic field. We see that as expected there are soft modes around ±kp = ^-j L . 
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FIG. 4: Boundaries of the holon-spinon excitation continuum in zero magnetic field for (a) it = 1 and density n c = 0.5 and (b) 
for u — 1 and n c = 0.2. 



B. Holon - 3 Spinon Excitation 

Other kinds of excitations with the quantum numbers (23), (22) involve more than two "elementary" excitations. 
As we are interested in the threshold of the spectral function, we focus on excitations that lead to the smallest possible 
energy for a given momentum. This leads us to consider holon - 3 spinon excitations characterized by a solution of 
the Bethe ansatz equations (7) with 

N = N GS - 1 , M = AI GS - 2. (30) 

It follows from (9), (10) that both Ij and J a are half-odd integers. In order to see that we are dealing with a 
four-parametric excitation we consider the number of vacancies for the integers Ij and J a (11). As 



-I < s f • 



there is the same number of vacancies for the Ij's as in the ground state, but we have one fewer integer, i.e. one 
additional "hole" . Similarly we have 

\J a \ < l(N GS -M GS ), (32) 

which tells us that we have one more vacancy than in the ground state but two roots less, which amounts to three 
"holes" . The energy and momentum of this excitation are 



E hsss = -e e (fc fc )-£;e a (A£) , 
i=i 

Phsss = -Pc(k h )- ^>( A t)j . ( 33 ) 

where the dressed energies and momenta are defined in (19) and (28), respectively. We note that there is no aditional 
constant contribution to the momentum because both the Ij and J a are half odd integers, leading to symmetric 
distributions of roots in the absence of holes. In Fig. 5 we plot the upper and lower boundaries of the holon - 3 spinon 
continuum for two different band fillings in zero magnetic field. We see that in part of the Brillouin zone the holon - 
3 spinon threshold is lower in energy than the holon-spinon continuum. 
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(a) (b) 

FIG. 5: Boundaries of the holon - 3 spinon excitation continuum in zero magnetic field for (a) it = 1 and density n c — 0.5 and 
(b) for u = 1 and n c = 0.2. The boundaries of the holon-spinon continua are depicted by the dashed blue lines. 



C. Thresholds in Zero Magnetic Field 

Wc are now in a position to determine the lower boundaries of the holon-spinon and holon - 3 spinon continua in 
the vicinity of kp. We find that the lower boundary is obtained as follows: 

• < P hs < k F 



Here the threshold is identical for both types of excitation considered above. The leading singularity of the 
spectral function is given by contribution of the holon-spinon continuum. The threshold is obtained as follows: 
we choose the plus sign in (27) and place the holon at its right "Fermi point" , i.e. 

I h = ^f- , k h = Q. (34) 

It therefore carries zero energy and momentum — nn c . For P/,. s = kp the spinon has rapidity A' 1 = — oo, which 
corresponds to momentum kp = ^rp. Taking the spinon rapidity from A h = — oo to A h = then traces out the 
lower boundary of the holon-spinon continuum for < P^s < kp. 

• 2k F > P hs > k F 

Here the threshold is given by the holon - 3 spinon continuum. It is obtained by placing the holon at its left 
"Fermi point" 

/" = -^§^, fc* = -Q. (35) 

It therefore carries zero energy and momentum 7rn c = 2k p. Two of the spinons arc placed at their left and right 
"Fermi points" respectively 

Jx,2 = , Af >a = ±oo. (36) 

Their contributions to the total momentum cancel. To obtain total momentum Ph sss = kp the third spinon is 
taken to have rapidity A' 1 = oo, which corresponds to momentum —kp = —^-Sr- Taking the spinon rapidity from 
A /l = oo to A h = then traces out the lower boundary of the holon - 3 spinon continuum for kp < Phsss < 2kp. 

At this point an obvious question to ask is whether other excitations may lead to even lower thresholds. As we 
have restricted ourselves to commensurate band fillings and zero magnetic field the answer close to kp is negative. 
This follows from the general expression for the low-energy part of the spectrum (5). This strongly suggests that the 
threshold of the spectral function is given by the two excitations given above and we will assume this to be the case 
in what follows. 
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D. Holon-Spinon Threshold in Zero Magnetic Field 

It is conceivable that the spectral function A(w, q) could exhibit features associated with the theshold of the holon- 
spinon excitation in the momentum range where the threshold of A(w, q) occurs at a lower frequency. We therefore 
determine the threshold of the holon-spinon excitation in zero field in the momentum range 2kp > Py ts > kp. It is 
obtained by choosing the plus sign in (27) and place the spinon at its "Fermi point" , i.e. 

Jh = _Ngs_2 Ah = _^ (3?) 

It therefore carries zero energy and momentum For P^ s = kp the holon has rapidity k h = Q, which corresponds 
to momentum —nn c . Reducing the holon rapidity from k h = Q traces out the lower boundary of the holon-spinon 
continuum for kp < Py is < 2k f. 



IV. THRESHOLD AT < P < k F FOR THE HOLON-SPINON EXCITATION 

In the following we consider the threshold at < P < kp for the holon-spinon excitation. By virtue of (34) we 
are dealing with a distribution of integers Ij that is symmetric around zero. The free parameter is then the spinon 
rapidity A h . More precisely, the distributions of integers are 

Ij = , l<j<N GS -l, 



Mgs + 1 



Q 



if 1 < a < ^f±i + J h 



J* - |_Mc|+i +a + 1 ii Mc i±1 + J h< a < MGs _ 1 - (38) 
Our goal is to determine the finite-size corrections to the energy of this excitation in the limit 

L,J\N GS ,M GS ^^ , L^^MSf* fed. (39) 



A. Finite-Size Corrections 



Our starting point are the Bethe ansatz equations (7) for the holon-spinon excitation where the holon sits at the 
Fermi momentum of the fcj's. It is convenient to write them in terms of counting functions 



z c (k) 



M 



sin k — A r 
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1 M 



A- A 
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The Bethe Ansatz equations then read 



z c{kj) 



2nJ n 



L 1 sv L 
We now turn these into integral equations by means of the Euler-Maclaurin sum formula 



1 ™2 

n— ri\ 



dx f(x) 
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24L 2 



where 



(40) 
(41) 

(42) 
(43) 



n + = n 2 + 



ri- = ni 



(44) 
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For the specific exitation we are considering we have 

T N GS -1 
1+ = ±- 



J± = ± 



2 2 

Using (43) in (41) and then changing variables from x (i.e. integers divided by L) to the rapidities we arrive at 

sin k - A \ 1 „ / sin k - A y 



(45) 



*(*) = k+ dA p.(A) - S"^) 
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.(A) 



24L 2 



ai(sinfc — ^4 + ) ai(sinfc — 



dfc 



Ps{A + ) 
A — sin k s 



o(L" 



/A-smfc\ .,, i A + JA , .... „/A-A'\ l./A-A£\ 



24L 2 



ai(A — sin<5 + ) cosQ + a x (A — sinQ_) cosQ- a 2 (A — A + ) a 2 (A — A_) 



Pc(Q+) 



Pc(Q-) 



(46) 



+ o(L- 2 ). (47) 



Here a n (x) is given in (17), the integration boundaries are fixed by 



Zc{Q±) = '—J^- i z s (A ± ) 



and the root densities arc defined by 



2np s (A) = , 2vp c (k) 

In addition there is the equation fixing the position of the hole 

=,(AJ) - 2 ' J " 



2ttJ ± 



dz c (k) 
dk ' 



(48) 



(49) 



(50) 



Here our notation makes the L-dcpendcncc of the rapidity of the hole explicit. Taking derivatives of (47) we obtain 
coupled linear integral equations for the root densities p cs 



2n 

cos k 



p c (k) = h cos(fc) / dA p s (A) ai(sin k — A) 



A- 



ai (sin k — A*l) cos k 
L 



24L 2 



a[ (sin k — A + ) a[ (sin k — A. 



Ps{A- 



p s (A) = / oi(A - smk)p c (k) 



dA' p s (A') a 2 (A - A') + -a 2 (A - A h L ) 

1j 



24L 2 



a[(A — sinQ + ) cos(J+ a [(A — sin Q_) cos Q_ a' 2 (A — ^4 + ) a' 2 (A — A_ 



Pc(Qh 



p«(^ 



P. (-4-) 



In shorthand notations these can be written as 



(51) 



(52) 



(53) 



where * denotes convolution on the interval [Q_,Q + ] ([A_,A + ]) if /3 = c (j3 = s). The components of the matrix 
kernel are 



K cc {k,k') = 0, K cs (k, A) = cos(fc) ai(sinfc — A) 

K sc (A,k) = ai(sinfc — A), K SS (A, A') = -a 2 (A - A'). 

It is useful to introduce a unified notation for the integration boundaries 



x± = 



Q± if a = c 
A± if a = s. 



(54) 



(55) 
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Using the Eulcr-MacLaurin sum formula on the expression for the energy (12) we obtain 

(0)' 



E = Lu + Lj2 dz e(°)(.) p a (z) - 4 0) K) 

^ Jx a 24i Pa 



" (X«) ej 0) '(Xg) 
(X?) p a (X«) 



(56) 



where the prime denotes the derivative with respect to the argument and ei are the bare energies, i.e. the driving 
terms in (19) 



4 0) (*0 = -2 cos(fc) — p — 2u — B , e^ ' (A) = 25. 



(57) 



Now we have to solve the system of equations (52) with boundary conditions (49) in by expanding both the densities 
and the integration boundaries in inverse powers of L. As the integral equations are linear we may proceed by first 
keeping the integration boundaries general and only expanding the densities as 



p a (z) = p a ,o(z) + jP a ,i{z) + 2^j2 



The various parts then fulfil the integral equations 



?a,a = Pi 0) a + K a p * pp. a , a = 0, 1 



/a/3 — ^% +K aJ * f lf3 



O) 



where the driving terms are given by 



2ir ' 

p<®(z) = -K as (z,A h L ), 



d 



dz 1 



K a p(z,z'). 



(58) 

(59) 
(60) 

(61) 
(62) 

(63) 



The integral equations for the /'s are solved by formal (matrix) Neumann-series and using this we can bring the 
expression for the energy in the form 



E = Lu + L^[ + dze^(z) Pa ,o(z) ~ $H*l) +Y, I + dz (z) p aA (z) - -^-(v s + v c ), (64) 
where the spin and charge velocities are 

(65) 



2irp afi (X a ) 



a = c. s. 



Let us denote the first two terms in (64) as LeGs({X±}). We consider it as a functional of the integration boundaries 
and expand it to second order around ±X a , i.e. 



e GS ({X£}) = e GS ({X a }) + Y, 

4e 



/9,<T 

s 2 



sx^ 



e GS ({XZ}) 



X%=<jXf< 



5X^5X7 



x%=„xf> 

X^=tX1 



ecs({A^}) 



(xP-<TxP)(xy-TXf). 



(66) 



We find that the linear term vanishes by virtue of the equations (19) for the dressed energies. For the quadratic term 
we find after some calculations 



eGs{{X±}) = S a {35<j T 2irv a [p a ,o{X a )Y 



(67) 
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The third and fourth terms in (64) can be simplified using the integral equations for the dessed energies. Introducing 
the shift of the hole rapidity in the finite volume by 



A h L = A h + ySA h , 



we can express the finite-size energy as 

E = Le GS ({X a }) - e s (A h ) + Lit ]T v a { [p a , (X a )(X™ - X a )] 2 + [p afi {X a )(X« + X a )] 2 } 



We may calculate SA h from the equation 
2irJ h 



z s (A h L , 



Q+ , s /A'/ -sinfc 

dk p c , {k) 9 1 



A + / \h _ A 

dA Ps , (A) e 1 



2 II 



1 



Q 



+ j- I dk PcA (k) 1 



A 1 } - sinfc\ 1 



u 



l\ dA PsA (A) 



2u 



Using the definition (68) we obtain 

L 



SA h 



2ir Ps , (A h ) 
L 

'2v Psfi {A h ) 



— / dk p Ci i(k) 



A h - sin k \ 1 



--/ dAp sA (A)6 



A h - A 
2u 



Here 



^;\A h ) 



app, Q {X^)g^{A h ) + C dk r%\k)6 f 



A dAr%\A)~ 

■A 



2u 



where g c<cx (A h ) = 8^- — "^WA g s rT (A h ) = — 9^ A '* 2n ' Tj4 ^ and the functions r^g fulfil the integral equations 



(68) 



(69) 



+ o(L- 1 ). (70) 



(71) 



(72) 



(73) 



B. Relating X" — aX a to Quantum Numbers 

In the next step we want to express the deviations of the integration boundaries from their thermodynamic values 
through appropriate quantum numbers. We define N c>s by 



N n I+. - I- 



L L 

N„ J + - J- 



L 



+ 

dk p c (k) , 
Q- 

dA Ps (A) . (74) 



We note that the number of down spins is A^ s — 1 rather than N s as we have one "deep" hole in the distribution of 
A's. The other quantities we want to use are 

2D C = 1+ + I_ = — [z c (Q+) + z c (Q_)} , 

2D S = J + + J_=J± [z s (A+) + z s {A_)\ . (75) 

Z7T 
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Using the integral equations (47) for the counting functions z c _ s we can rewrite these as 

"A- />oo 

2d, = / dA Ps (A) - / dAp s (A), 



2D S 
L 



2D C 
L 

We have to order 1/L 



2d, 



dk p c {k) — I dk p c (k) — 



1 f A + 




- dA0 


ft 


I" J A- 





ft (A) 



-A- 



x + 1 r x " 

dzp afi (z) + — dz p a .i(z). 



x° 



The second term no longer depends on X£ and is denoted by 



A, 1 



dk p c ,x{k) , 



dk p s ,i(A) 



-A 



The variation of the integration boundaries X" with np is now easily calculated to leading order in L 



dXl 1 Z, 



a/3 



Here the dressed charge matrix Z 



Z 



drip 2p a .o(X a ) 



is given by Z a p = £ a p(XP), where £ a p fulfil the set of coupled integral equations 



Similarly we have 



2d, = 



r°° 2 ■ 

dA Ps , (A) - / dA Psfi (A) + -£>™p, 

i JA + L 

fQ- fir I [A+ / fo\ 

2d c = / dk p Ct0 {k) - / dk p cfl (k) - - / dA - p s , (A) + —D 

J-n Jq + kJ a _ \uj L 



imp 

c I 



where 



2D uap 



dA/9 Si i(A) - / dAp, M (A) , 



2D 



imp 



-Q 



dk p c ,i(k) 



1 f A /A\ 1 fA h 

dk p cA (k) - dA 9 - p s ,i(A) + -0 — 



(76) 



(77) 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



These allows us calculate the dependence of the integration boundaries on d a 

dX± Zp a 
ddp ~ Pa,o{ Xa )' 

Combining (79) with (84) we obtain the described relationship between the change in integration boundaries and the 
quantum numbers n a , d a 



(84) 



X^TX a = ± 



1 z, 



2p afi {X a ) 



1 

1"* 



An - ^N™ p 



7 T 



P a ,o( Xa ) 



Ad, - ±D»p 



(85) 



Here An C)S and Ad c s are the changes in the quantum numbers (74), (76) compared to the ground state. 
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C. Result for the Finite Size Energy 

Putting everything together we then arrive at 
E = Le GS ({X a }) - e s (A h ) 

e' s (A h )SA h } 



1(1"/ s „ 

Z I ~6 (Ws +Vc) + 



\AN 7 (Z T )-i Va Z-iANp + D 7 Z ia v a Zl p D p 



Here we have defined 



AiV Q = AN a - N™v , 
D a = D a - D™p . (87) 

In writing (87) we have used that for the ground state D a = 0. The first term in (86) is the ground state energy, the 
second term the excitation energy in the thermodynamic limit and the second line gives the L _1 corrections. For our 
case we have 

N c = 1+ - I- = N GS - 1 , N s = J + - J_ = M GS , D c = = , D s = J+ ± — = . 

For the ground state we have 



N c = N GS , N s = M GS , (89) 



which gives 



AN C = -1 , AN S = 0. (90) 

D. Relation of iV), mp to the spectrum 

The quantities 7V™ P are given in terms of the solutions to the coupled integral equations (59). It is possible to 
relate them to properties of the dispersions of the elementary excitations as follows. The integral equations for p a .\ 
are formally solved by 

p„,i - (i - k)i\ * P f x = -(i - k)~l * k Ps . (91) 

Substituting this into the equations for 7V™ P we obtain 

NlT P = — ( Q dk (1 — K)-J * K P „ = - [ Q dk (1 - K)-\k, A h ) , 
J-Q J-Q 



NT P = - J dA (1 - K)-£ * Kp s = 1 - J dA (1 - K)- 1 (A, A"). (92) 

Here we have used e.g. 

(i - k)-J * k Ps = (i - kyj *(k-i + i) 0s = (i - k)- 1 . (93) 

On the other hand, by differentiating the integral equations (19) for the dressed energies we obtain 

de a dtp y 

= ^ * K, a - S ac , (94) 



which are solved by 



This gives us our first relation 



%r = - f Q dk(i-k)-^(k,z a ). (95) 

OH J-Q 



jyimp = d6s ^ H \ (96) 
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By comparing (94) to (81) we observe that 



de a 



(97) 



which in conjunction with (96) allows us to relate iV™ p to the dressed charge matrix as iV™? = -£cs(A h )- This 
generalizes the analogous relation for the spin-1/2 XXZ chain found in Ref.[18]. A second relation is obtained by 
considering 



dB ~ dB 



* Kfi a + 25 as — 5 ac 



Comparing this to (81) and using the linearity of the integral equation we observe that 

de a 



dB 



2£sa £e 



The formal solution of (98) is 



Our second relation is then 



de a _ cka 
dB ~ d f i 



+ 2 dk(l-K)- s l 

J -A 



N imp = l __ 



de s (A h ) 8e s (A h ) 



dB 



Of i 



= i-6,(A /l ). 



(98) 



(99) 



(100) 



(101) 



E. Simplification for zero Magnetic Field 

In the absence of a magnetic field we have A = oo which allows us to simplify our results for the finite-size corrections 
to the energy (86). The dressed charge matrix takes the form 1 

z= a °i ) , (102) 



,2 ^2/ 

where £ = £(Q) is obtained from the solution of the integral equation 

rQ 



£(fc) = 1 + / dk' cos(fc') i?(sin(fc) - sin(fc')) £,(k'). (103) 

J-Q 

Here the function R(x) is 

r°° duj e lu)x 

R W = / ^ TT ro i iv ( 104 ) 

J_ 00 Z7r 1 + exp(2u\uj\) 

and can be expressed in terms of the Digamma function. The integral equations for the dressed energies and root 
densities simplify to 1 

f Q 

e c (k) = -2cos(fc) - n - 2u+ / dk' cos(fc') R(sm(k) - sin(fc')) e c (k') , 

J-Q 

rQ 

e s (A) = / dkcos(k) s(A - sin(fc)) e c (k) , 



1 f Q 

p c (k) = hcos(fc) / dk' R{shx(k) - sin(fc')) Pc(k ) , 

2tt J_ q 

rQ 

p s {A) = / dk s(A - sin(fc)) p c (k) , (105) 
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where 



s{x) 



Au cosh ( 5^ 

2u 



The finite-size energy can be expressed as 



E = Le GS ({X a })-e s (A h )-^-(v s +v c )-ye' s (A h )8A h 



2ltV c 



(AN C - iV™?) 2 

if 5 



2ltV s 



2 I 



51 
2 



where 



,yimp 

rjimp 
^ s 

2D'T P 



2iv imp _ 1 
. 



-0 



dkp Cjl (k) 



dk [p c ,i(-fc) - Pc,i (A;)] + ^0 (tanh 

r(H^)r(i-»f) 



dk p c ,i(k) i In 



; sin k 1 



Here p c i is the solution to the integral equation 

cos(/c) 



Pc,i(k) 



4 M cosh r (Ah - Binfc) 



cos(fc) / dk' R(sink — sink') p Ci i(k'). 



(106) 



(107) 

(108) 
(109) 

(110) 
(111) 



The fact that D" np = is established in Appendix A. The relation (96) of TV™ 13 to the dressed energy of spin 
excitations remains valid for zero magnetic field, i.e. we have that A™ p = d£s g^ ^ ■ The quantum numbers AN a , D a 
for the holon-spinon excitation in zero magnetic field are 



D c = D s = , AN C = -1 



AN, = 0. 



(112) 



This shows that there are no L _1 corrections to the excitation energy in the spin sector. In terms of the field theory 
picture (see section VI)of a deep hole interacting with the gaplcss Luttingcr liquid degrees of freedom this shows that 
the interaction of the deep hole with the gapless spin sector is indeed irrelevant. Due to the presence of the marginally 
irrelevant interaction of spin currents we expect logarithmic corrections, but these are beyond the accuracy of our 
finite-size calculation. 

Explicit values for the quantities N^, mp and D™ p are readily obtained from a numerical solution of the linear integral 
equation (111). We present results for several values of u and two band fillings in Figs6 and 7. 



1. Large-u Limit 

In the strong interaction limit u 3> 1 we can solve the integral equation (111) by iteration and obtain explicit 
results. To leading order we find 

g im P _ Ps{A h ) _ P hs N - wv _ sin(Q) 

c ~ 2im c 4fc F ' c 2ucosh(7rAV2u)' 



V. THRESHOLD AT k F < P < 2k F FOR THE HOLON - 3 SPINON EXCITATION 



We now want to determine the finite-size corrections to the energy of the holon - 3 spinon excitation along the 
low-energy threshold in the momentum range kp < P < 2kp. The threshold of this excitation is described by the 
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(a) 



-0.2 
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0.25 



0.5 

P/K 
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u=1 

u=2 



0.75 





-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 



(b) 



0.25 



0.5 

P/K 



u=0.5 

u=1 

u=2 



0.75 



FIG. 6: A^ mp as a function of momentum for u = 0.5, 1, 2 and densities (a) n c = 0.2 and (b) n c = 0.5. 





FIG. 7: DJ. mp as a function of momentum for u = 0.5, 1, 2 and densities (a) n c = 0.2 and (b) n c = 0.5. 



Bcthc Ansatz equations (7) with N = N G s — 1, M = M G s — 2 and half-odd integer numbers 



iV GS - 1 

Ij = ^ + 3 , 1 < j < N GS - 1 



2 



M G S i Th 
J h <a< Mi 



if 1 < a < + J 



GS 



(114) 



The corresponding distributions of half-odd integers are shown in Fig. 8. The calculation on the finite-size corrections 
to the energy now proceeds just like for the holon-spinon excitation. The result is given by (86), (87), (78), (83) where 
now 

1 



Nr. 



/_ = N GS - 1 , N s = J+ - J_ = M, 



GS 



1 , Dr. 



D K = 



J+ + J- 







Recalling that in the ground state the N c ^ s "quantum numbers" arc 

N c = N GS , N s = M GS 
we conclude that for the holon - 3 spinon excitation we have 



-1 . AN, 



1 , D n 



1 

2 ' 



D s =0. 



(115) 
(116) 
(117) 
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J h 



FIG. 8: Distribution of integers for the threshold of the holon - 3 spinon excitation. The hole in the distribution of Ij is at 
— (Ngs — l)/2 and the three holes in the distribution of J a occur at ±Mgs/2 and J respectively. 

The reduction in the zero magnetic field case is completely analogous to the holon-spinon excitation. Hence the energy 
is again given by (107), (110), (111), but the "quantum numbers" N c>s , D c s are now given by (117). 

VI. FIELD THEORY 

We now relate our results to the field theory treatment of Ref. 22 for the threshold singularity problem. There 
it was shown that a high-energy excitation in a spinful Luttinger liquid can be mapped to a mobile impurity in a 
Luttingcr liquid. The corresponding Hamiltonian is H = ^2 a=c s H a + Hd + Hi nt , where 



H a = — I dx 



1 fd<S> a 



K a 



fdG a 



K a \ dx J \ dx 

H d = I dx dfi(x) [e s {P) - iud x ] d(x) , 



H 



int 



= I dx 



V R -V l V r +Vl - 
o x Q c d x <S> c 



2tt 



2tt 



d(x)d\x). 



(118) 



Here the Bose fields $ Q and the dual fields Q fulfil the commutation relations (4), d{x) and cr are annihilation and 
creation operators of the mobile impurity, which carries momentum P and travels at velocity u. The parameters Vr,l 
characterize the interaction of the impurity with the low energy charge degrees of freedom. The parameters of H c s 
and Hd in (118) arc readily identified with quantities obtained from the Bcthe Ansatz. The spin and charge velocities 
v SfC are given by (65) and the Luttinger parameters are 



Ks = 1 , K c 



(119) 



where £ is given by (103). The velocity of the impurity is expressed in terms of the solutions to the integral equations 
(105) as 



e>(A h ) 



(120) 



2n Ps {A h )- 

The "chemical potential" of the impurity is e s (A h ), where the position A h of the hole is fixed by the requirement 

P hs (^ h ) = P (121) 

The parameters Vr,l entering the expression for H lnt are deterined as follows. Following Refs 18,22 we remove the 
interaction term H lnii through a unitary transformation on the fields 



U = cxp I —i I dx 



2tt 



(122) 



where 



-(V l tV r )K c 



(v c + u)A(5_, c ± (v c — u)AS 



+,c- 



(123) 
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In the resulting Hamiltonian the impurity no longer interacts explicitly with the charge part of Luttingcr liquid, but 
it does affect the boundary conditions of the charge boson. In particular we find that 



d x $ c 



u^d x ^ c u = d x $ c 



d x e c = u^d x e c u = d x e c 



iKc. 



1 



AS 



A6_, C W , 



AS 



AS- 



2VK C 

Taking the expectation value of (124) with respect to a state with a high-energy hole we find that 

1 



(124) 



dx (d x Q c ) = J dx (d x Q c 
dx (d x $ c ) = dx (d x $ c 



2VK, 

s/kz 



A<5+. c + A<5_, c , 
A5 +c - AS- 



(125) 



Denoting by p Ci R and p c ^ the charge densities at the right and left Fermi wave numbers respectively we have that 
the numbers AN and D of low-energy charge and current excitations is 



AN — J dxV2 [p c>R + p c>L ] = -- -H j dx (d x $ c ) , 
AD = s/2 j dx [p c , R - Pc . L ] = ^ J dx (d x Q c ) . 



(126) 



We are now in a position to identify the additional contributions in (125) that arise from the interaction of the 
Luttingcr liquid with the mobile impurity with the quantities iV™ p and £>™ p . 



A. Holon-Spinon Exitation 

Taking into acount that for the holon-spinon excitation we have 

-— [dx (d x $ c ) = -1 , — [ dx (d x Q c ) = -1 , 

T J 7T J 

we conclude that the quantities iV™ p and 2Z?J, mp are related to the phases AS± yC by 

A<L_ r - A5_ r 



(127) 



N imp 



2D\ 



■y/2K, 



1 



2tt 

1 AS +tC + AS- 



\f2K c 



2 71 



(128) 



In Ref. 22 it was shown that in the momentum range \P\ < the single particle spectral function exhibits a threshold 
singularity of the form 



A(uj,P)oc(uj-e s (P))-^- 
where the exponent is expressed in terms of the phase-shifts A5± jC by 



(129) 



Mo, 



IKc , A£ +)C + A£_, 
2 



2tt 



V2K C 



AS+. C - A<5_, 
2^ 



= 1 - 



N 1 - 



We note that (130) differs from the Luttinger liquid result 6,23 

i Kc 

1-lLL = 1 T ' 



1 

AK< 



(130) 



(131) 



However, in the zero energy limit P — > kp the Luttinger liquid result is recovered. Explicit results for p,o - are 
obtained by solving the integral equation (111) numerically. In Fig. 9 we present results for density n c = 0.6 and three 
different values of u. 



19 




FIG. 9: Threshold singularity exponent ^o,- as a function of momentum for n c = 0.6. 



1. Large-u Limit 



For u> 1 wc have K c h and using (113) we find 



Mo.- 



1 1 

2 ~ 8 



P 



(132) 



This suggests that the singularity is only weakly momentum dependent and is close to a square root. The result 
(132) agrees with the exponent computed by exploiting the factorization of the wave function into spin and charge 
parts 28-30 . 

2. Comparison to DMRG Results 

In Ref.4 the single particle spectral function for parameter values u = 1.225 and n c = 0.6 was computed by the 
dynamical density matrix renormalization group method 5 . The following two exponents were reported based on a 
scaling analysis of the low-energy peak heights 



Mo, 



0.86 for P = , 



Mo, 



0.78 for P = —. 

10 



(133) 



The Luttingcr liquid parameter and Luttinger liquid threshold exponent are K c = 0.6851 and Ull = 0.4638 respec- 
tively. Our results are 



Mo,- 



0.864 for P = , 



Mo, 



0.832 for P = 



10' 



(134) 



The agreement for P = is excellent, but the P = jq values are very different. One possible explanation for this 
discrepancy is as follows. The frequency range over which the power-law (130) holds can a priori be quite small and 
furthermore is expected to narrow as P increases from zero to kp. Extracting the threshold exponent from the scaling 
of the peak height could then require very small values of the parameter r\ in Rcf. 4. 



B. Holon - 3 Spinon Exitation 



The threshold of the holon - 3 spinon excitation can also be analyzed in terms of a mobile impurity field theory 
model 12,22 . The threshold behaviour in the range kp < P < 3kp is no longer singular but describes a power law 
"shoulder". The exponent can be expressed in terms of the quantities D" np , 7V), mp derived above. 
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VII. ELECTRON GAS 



The analysis for the electron gas is completely analogous to the one for the Hubbard model. The Hamiltonian is 



N 



i:) 2 



H = - ]T — - + 4uJ2 S(xi - x 3 ) ~fiN + B(2M - N). 

3=1 i<3 



(135) 



The Bcthc Ansatz equations read 



26,27 



A I 



kjL = 2nlj — 



kj A a 



, j = l,...,N 



2 ii 



(136) 
(137) 



For the ground state and the excitations we are interested in we only need to consider real roots of (136), (137). For 
these the "quantum numbers" Ij , J a fulfil the restrictions 



Ia is 



integer 



if M is even 



half — odd integer if M is odd, 



(138) 



J a is 



integer 



if N - M is odd 



half — odd integer if N — M is even, 
The energy and momentum of such Bethe ansatz states are 



N 



N 



E = 2BM + J2i k 



3=1 



B], P = J2 k 3 

3=1 



2tt 

T 



N M 
3=1 



o=l 



(139) 



(140) 



A. Ground State and Excitations in the thermodynamic limit 



In the thermodynamic limit the ground state is described by the solution of the integral equations 

1 f A 

Pco(*) = 7^-+ J A dK Kcsik- A) p S:0 (A) , 



P«,o(A) 



dk K SC (A - k) Pc . (k) + / dA' K SS (A - A') p, i0 (A') , 

-Q J -A 



where we have defined integral kernels as 



K cc {k,k') = 0, K cs (k,A) = ai {k-A) , 

K sc (A,k) = oi(fe-A), AC„(A J A / ) = -a 2 (A-A'). 



Density and magnetization per site are given by n c and % — n s respectively, where 

rQ 
'-Q 



rQ rA 
/ dk p c fi(k), n s = / dA p a ,o(A). 
J-q ^ —A 



(141) 



(142) 



(143) 



The dressed energies of elementary charge and spin excitations in the thermodynamic limit at finite density and 
magnetization are given in terms of the solutions to the coupled integral equations 



fA 

e c (k) = k 2 -p-B+ / dA )C cs (k - A) e s (A) 
J -A 



rQ fA 
e,(A) = 2B+ I dk IC SC (A - k) e c (k) + / dA' IC SS {A - A') e s (A') , 

Q J -A 



(144) 
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where the integration boundaries are fixed by the requirements 

e c (±Q) = , e s (±A) = 0. (145) 

The holon-spinon excitation is constructed in complete analogy with the Hubbard model. In the thermodynamic limit 
its energy and momentum are expressed as 

E hs = -e c (k h ) - e s (A h ) , P hs = - Pc (k h ) - Ps (A h ) ± nn c . (146) 



B. Finite-Size Corrections for Threshold Exitations at |P/i S | < kF 

Following through the same steps as for the Hubbard model we obtain 
E = Le GS ({X a }) - e s (A h ) 



7T 27T 

6L {Vs+Vc) + T 



- ye' s {A h )SA h + oiL- 1 ) . (147) 



Here eGs({^ Q }) is the ground state energy per site in the thermodynamic limit (and we again use notations where 
X c = Q and X s = A), —e s (A h ) is the 0(1) contribution of the spinon excitation in the thermodynamic limit and 
v s and v c are the velocities of the gapless elementary spin and charge excitations. They are given in terms of the 
solutions of the integral equations (141), (144) by 



27rp c ,o(<3) ' ^Ps,o{A)' 
The dressed charge matrix Z in (147) is defined as 

tcc(Q) ics{A) 



(148) 



z = [mm • (149) 



where fulfil the set of coupled integral equations 

rX 
-X 



/X^ 
dz 1 £ ay (z 7 ) /C 7/3 (z 7 - zp) . (150) 



The quantities A7V Q and D a are defined as in the case of the Hubbard model 

AN a = AN a - N™v , D a = D a - D™? , (151) 

where now 

AT P - [ Q dk PcA (k) , N™v = [ A dA p sA (A) , (152) 

J-Q J-A 

/—Q />oc /*— A poo 

- dk PcA (k), 2£>f p = / dAp.,i(A). (153) 

-oo JQ J — oo J A 

Note that these are different from what we had for the Hubbard model. The root densities p C; i and p Si i fulfil the 
coupled integral equations 

Pci(fc) - - ai (k-A h ) + / dA/C„(*-A)p a ,i(A) , (154) 

J-A 

Ps.i(A) = a 2 (A-A fc ) + / dkK sc {A-k) Pc .i{k) + / dA' fC ss (A — A') p a ,i(A') . (155) 

J-Q J-A 

Finally, e^,(A) is the derivative of the dressed energy (144) and SA h /L is the shift in the hole rapidity due to the finite 
volume quantization conditions. It is obtained from a set of equations completely analogous to (71), (72) and (73). 
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C. Simplification for zero Magnetic Field 

In the absence of a magnetic field we have A = oo, which again allows us to simplify all expressions. The integral 
equations for the dressed energies can be written in the form 



e c (k) = fc 2 -/i+ / dk' R(k' - k) e c (k') , 



-Q 



*.(A) 



dk s(A - k) e c (k) , 



(156) 



where the integration boundary Q is fixed as a function of the chemical potential /i by the requirement e c (±Q) = 0. 
The dressed charge matrix takes the form (102), where £ = is obtained from the integral equation 



=1+1 dk' R{k - k') £(£;')• 



(157) 



The expression for the finite-size energy simplifies to 

E = Le GS ({X a }) - e s (A h ) - ^.(v s + v c ) - ±e' s (A h )8A h 



2nv c 
L 



(AN C - iV™?) 2 
4^2 



D, 



27TU, 



(an s - ±an c - 



Pi 
2 



o(L" 1 ). (158) 



Q 



The expressions for A^ mp and D™ p become 

iV™P = 27V] m P-l= / dk PcS {k) , 

J-Q 

2D imp = rdklpdi-Q-pdik)], D\ m ? = 

JQ 

where p c .i is the solution to the integral equation 

1 



Pc,i(k) 



Au cosh ( ^ 



+ / dk' R(k-k') p c ,i{k'). 



(159) 
(160) 

(161) 



Here R(x) is given by (104). We note that our expression for iV™ p is related to the dressed energies (144) by 

de s (A h ) 



N nnp 



dp 



(162) 



For the holon-spinon excitation we have 

AN C = 
AN S = N s - 



N c - N GS = (N GS - 1) - N GS 



(163) 



Like fo the Hubbard model the corrections in the low-energy spin sector vanish, which in terms of the field theory 
description implies that the interaction of the high-energy spinon with the low-energy spin sector is irrelevant. 



D. Threshold Exponent 

As for the Hubbard model the spectral function exhibits a threshold singularity of the form 

A(uj,P)<x(uj-e s (P))-^- , 



(164) 
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where the exponent /^o,- is expressed both in terms of the phase-shifts S±, c arising in the field theory treatment (122) 
and in terms of the quantities JV™p, L»™ p in (160) 



1 

l -2 




J +,c 



-i 2 



2- 



1 



S-.c 



2- 



iV lmp 



(165) 



In Rcf.22 the phase-shifts (5± jC were expressed in terms of properties of the excitation spectrum. We have verified the 
relations given there by numerically solving the relevant integral equations for the Yang-Gaudin model. 



VIII. SUMMARY AND CONCLUSIONS 



In this work we have determined the 0(L _1 ) corrections to energies of excited states in the Hubbard and Yang- 
Gaudin models for the case where in addition to any finite number of low-energy excitations, high-energy excitations 
are present as well. This extends the work of Woynarovich 23 , which considered exclusively low-lying excited states. 
We have focussed on the case of a single high-energy holon or spinon, but the method is easily extended to other 
cases. There are several contributions to the 0(L~ l ) energy corrections. One of these arises from the quantization 
of the momentum for the high-energy particle in the finite volume. More interestingly, we find that the presence of 
a high-energy particle leads to a modification of the low-energy part of the spectrum. This effect can be understood 
by a mapping to a model of a mobile impurity coupled to the spin-charge separated Luttinger liquid that describes 
the low-energy degrees of freedom 22 . By matching the finite-size spectrum of the mobile impurity model to the one 
obtained from the exact solution, we have obtained explicit results for threshold singularities in the spectral function 
A(u, q). In the momentum range \P\ < kp the negative frequency part exhibits a singularity at uj = e s (P) of the form 

/ \ -Mo,- 

A(u,P) cx \ uj-e s (P)j , (166) 

while the positive frequency part vanishes in a characteristic power-law fashion above uj = — e s (P) as 

A(w,P) oc (w + esCP)) 1 . (167) 

Expression (130) for /io,- is the main result of this work. The resulting spectral function is depicted in Fig. 2. 

There have been a number of previous studies of the single-particle spectral function of the one dimensional Hubbard 
model below half-filling. Ref. 4 reports dynamical density matrix renormalization group results for density n c = 0.6 
and U/t — 4.9. We found that the threshold exponent (130) is only in partial agreement with Ref. 4. A possible 
explanation is that smaller values for the imaginary part of the energy are required in the DMRG computation in 
order to extract the singularity reliably. It would be interesting to test this conjecture. In the large-it limit the 
single-particle spectral function has been computed 28 by exploiting the factorization of the wave function into spin 
and charge parts 29 . The behaviour obtained by this method was reported in Ref. 30 and agrees with (132). Finally, 
singularity exponents obtained by completely different methods have been reported in Ref. 30. We have checked 
that the numerical results for the exponent ^to,- in the range \P\ < kp, density n c = 0.59 and several values of u 
(u = 0.25, 1.225, 2.5) are in agreement with ours. It would be interesting to demonstrate the equivalence of the 
expressions for the exponents of Ref. 30 and our results analytically. 
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Appendix A: Zero Field Limit 

In this Appendix we consider the zero magnetic field limit for the quantity D™ p (83). By definition we have 

/>oo />oo 

D'r = / dA [p a ,i(-A) - p 8>1 (A)] = - / dA /(A), (Al) 

J A JO 
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where 



f(A)=p s , 1 (A + A)-p s , 1 (-A-A). 



(A2) 



After Fourier transforming the integral equation for p s> i we obtain the following set of equations for /(A) and p Ct -(k) = 
Pc,i(k) - Pc,i{-k) 



rQ 

/(A) = -R(A + A + A h ) + R(A + A - A h ) + / dk s(A + A- sink) p c ,-(k) 

J-Q 

/>OG 

+ / dA' [R(A - A') - R(A + A' + 2A)} f(A') , 



PcAk) 



Jo 

rQ 

cos(fc) [s(A - sinfc) - s(A h + sinfc)] + cos(fc) / dk' i?(sinfc - sinfc') p c -(k') 

J-Q 

/•OO 

cos(fc) / dA [s(A + A - sinfc) - s(A + A + sinfc)] /(A). 
Jo 



We now observe that for large A and |A''| <C A + A the driving term in (A3) is small 



R(A + A + A' 1 ) - R(A + A- A' 1 ) 



2u A h 



7T (A + A) 3 ' 

Iterating the integral equations in this limit then shows that Z?™ p vanishes when A tends to oo. 



(A3) 



(A4) 



(A5) 



Appendix B: Finite-Size Corrections for a high-energy holon excitation 

Our starting point are the Bethe ansatz equations (7) for the holon-spinon excitation where the spinon sits at the 
Fermi momentum of the A Q 's. The hole in the distribution of fcj's is denoted by k h and the correspoinding integer 
in the logarithmic form of the Bethe Ansatz equations (7) by I h . Expressing the Bethe Ansatz equations in terms of 
counting functions (41) we have 



Z c (fcj 



2nL 



where the integers Ij and J a are given by 



Jo 



Mi 



1 



2 2 

2 



Q 



Z S (Aq) 



a = 1, 



2ix,J a 



(Bl) 



GS 



if i <j < &m + i h 



(B2) 



Finite-Size Corrections 



As before we turn these into integral equations by means of the Euler-Maclaurin sum formula (43). This results in 



.(A) = 
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J A 
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rQ+ 
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sin fc — A 



1 



24L 2 



a\ (sinfc — A + ) a\(sva.k — A- 



Ps(A+) 



A — sin fc 



dA Ps {A') 9 



Ps(A-) 

(A- A' 



o(L~ 



V 2u 



Pc(Q+) 



Pc{Q-) 



Ps(A^ 



Ps(A-) 



(B3) 



o(L- 2 ), (B4) 



where a n (x) is given in (17), the root densities p CjS are given in terms of the counting functions by (49), and the 
integration boundaries are fixed by 



2irl± 



Zs(A±) 



2irJ± 



(B5) 
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Here 

N GS 1 
2 2 

The equation fixing the position of the hole is 



± 



r _ M GS 

— ^ ! J- 



9irJ~ h 

Zc(k h L ) = — = fixed. 



(B6) 



(B7) 



Here our notation makes the L-dependence of the rapidity of the hole explicit. Following through the same steps as 
in IV we find that the finite-size energy is expressed as 



E = Le GS ({X a }) - e c (k h ) 
If*"/ . 
+ L I ~6^ s +Vc ^ + 



-AN 7 (Z T )^v a ZjAN p + D, ( Z ia v a Z^Dp 



e' c {k h )5k h \, (B8) 



where eGs({A^"}) an d —e c (k h ) are respectively the ground state energy per site and the dressed energy of the holon 
in the thermodynamic limit, Z a p are the elements of the dressed charge matrix (80) and AN a , AD a are defined by 
(87), where 









primp _ 


/ dk p c ,i(k) , 


primp _ 
1 ¥ s 




J-Q 








poo 


2D hap = 
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/ dAp St 




J — DO 


J A 




r-Q 




2D hap = 


/ dk p Ci i(k) - 


/ dk p c . 




J — 7T 


Jq 



dk p s ,i(A) , 



-.1 



dAe(-) AM (A). 
Here the root densities p a> i fulfil the coupled integral equations 

Pd(k) = cos(fc) / dA ai(A - sin(fc)) p s ,i{A) , 
J -A 

rQ rA 
Ps,i(k) = -ai(A - sin(A:' 1 )) + / dk ai(A - sinfc) p c ,i(k) - / dA' a 2 (A - A') /3 s ,i(A') 

J-Q J-A 



The quantum numbers D a , AN a arc 



A7V C = 0, AN s = -l. 



Like for the spinon threshold it is possible to express -/V™ p in terms of the dressed energies. We find that 

„ , , de c (k h ) Nimp _ 1 \de c (k h ) de c {k h Y 



N" 



dp, 



dB 



dp 



(B9) 



(BIO) 



(Bll) 



(B12) 



2. Simplification for zero Magnetic Field 



In zero magnetic field the expression for the energy simplifies to 

E = Le GS ({X a }) - e c (k h ) - £-(v. + v c ) - ye' c (k h )Sk h 

6L L 



where 



2"KV r 



{AN C - N c 1TLp ) 2 
4£2 
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( AA r s -^ + \Y Dl 
2 2 



2D imp = 



Q 



i 

t 
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Q n J-Q 



: sin k ^ 



r(|-^)r(i + itr) 



(B13) 

(B14) 
(B15) 
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Here pc.i is the solution to the integral equation 



Pc,i{k) = - cos(fc)i?(sin(fc) - sin(fc' 1 )) + cos(fc) / dk' R{smk - sink') p c ,i(k')- (B16) 

J-Q 
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